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ABSTRACT 

For a robust interpretation of upcoming observations from PLANCK and 
LHC experiments it is imperative to understand how the inflationary dynamics 
of a non-minimally coupled Higgs scalar field with gravity may affect the deter- 
mination of the inflationary observables. We make a full proper analysis of the 
WMAP7+SN+BAO dataset in the context of the non-minimally coupled Higgs 
inflation field with gravity. 

For the central value of the top quark pole mass = 171.3GeV, the fit of 
the inflation model with non-minimally coupled Higgs scalar field leads to the 
Higgs boson mass in range 143.7 GeV < ran < 167GeV (95% CL) 
We show that the inflation driven by a non-minimally coupled scalar field to the 
Einstein gravity leads to significant constraints on the scalar spectral index ns 
and tensor-to-scalar ratio R when compared with a tensor with similar constraints 
to form the standard inflation with a minimally coupled scalar field. 
We also show that an accurate reconstruction of the Higgs potential in terms of 
inflationary observables requires an improved accuracy of other parameters of the 
Standard Model of particle physics such as the top quark mass and the effective 
QCD coupling constant. 

Subject headings: cosmology: cosmic microwave background, cosmological pa- 
rameters, early universe, inflation, observations 



1. Introduction 

The primary goal of particle cosmology is to obtain a concordant description of the 
early evolution of the universe, establishing a testable link between cosmology and particle 
physics, consistent with both unified field theory and astrophysical and cosmological mea- 
surements. On the ground, the Large Hadron Collider (LHC) at CERN is investigating the 
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elementary particle collisions in the TeV energy range, seeking to validate a large number of 
theoretical predictions of the Standard Model (SM) of particle physics and beyond. In the 
sky, the PLANCK Surveyor is actively taking precise measurements of the Cosmic Microwave 
Background (CMB) temperature and polarization anisotropies. 

Inflation is the most simple and robust theory capable of explaining astrophysical and 
cosmolo gical observations ; at the same time provi d ing self-consist e nt primordi a l initial con- 
ditions Jstarobinskv 1979 : buth 198ll : ISato 198ll : klbercht 19821 : iLinde 19821 iLinde 1983h 
and mechanisms for the qu antum generation of scalar (curvature) and tensor (gravita- 



tional waves) perturbation s ( IMukanov 198 It lHawking 19821 ; iGuth 1982c IStarobinsky 1982 



Bardeen 19831 ; lAbbot 1984 ). In the simplest class of inflationary models, inflation is driven 
by a single scalar field <ft ( or inflaton) with some potential V(4>) minimally coupled to the 
Einstein gravity. The perturbations are predicted to be adiabatic, nearly scale-invariant and 
Gaussian distributed, resulting in an effectively flat universe. 



The WMAP cosmi c microwave background (CMB) measurements alon e (IDunklev et al. 2 009 



Larson et al. 20101) o r complemented with other cosmological datasets (IKomatsu et al. 20091 ; 
Komatsu et al. 20101 ) support the standard inflationary predictions of a nearly flat universe 
with adiabatic initial density perturbations. In particular, the detected a nti-correlations be- 
tween temperature and E-mode polarization anisotropy on degree scales (INolta et al. 20091 ) 
provide strong evidence for correlation on length scales beyond the Hubble radius. 

Alternatively, one can look to the inflationary dynamics based on models beyond the 
Standard Model (SM) of particle p hysics. The hybrid inflation models involving supersym- 
metric (SUSY) TeV e nergy scales (IDvali et al. 19941 ) and minimal supergravity (SUGRA) 
( 
( 



Linde fc Riotto 19971) prov i de natural connection between cosmology and particle physics 



Cervantes-Cota fc Dehnenl ; ISenoguz fc Shafi 20051 ). The realization of these inflationary 



scenarios introduces new physics between the electroweak energy scale and the Planck scale, 



leading to distinct predictions of the main inflationary par ameters, such as the spectra 
n,s of scalar perturbations and the te nsor-to-scalar ratio R (Rehman et. al 20081; 12 009 



index 



2010). 



However, a number of recent papers (IBezrukcy fc Shaposhnikov 2 008; B arvinskv et al. 2 008 



Bezrukov et al. 2009: Bezrukov et al. 2009: De Simone et al. 2009 



Bezrukov et al. 20091 ) 



reported the possibility that the SM of particle physics with an additional non-minimally 
coupled term of the Higgs field to the gravitational Ricci scalar can give rise to inflation 
without the need for additional degrees of freedom to the SM. This scenario is based on 
the observation that the problem of the very small value of Higgs quadratic coupling re- 



quired by the CMB anisotropy data can 


:>e solved if the Higgs inflaton has a large coupling 


with gravity ( 


Futamase & Maeda 19891; 


Fakir & Unruh 199(1 


Komatsu & Futamase 1999; 


Tsujikawa k, Gumjudpai 2004; 


Barvinskv & Kamenshchik 1994 


)■ 



The resultant Higgs inflaton effective potential in the inflationary domain is effectively flat 
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and can result in a successful inflation for values of the non-minimal coupling constant 
£ ~ 10 3 — 10 4 , allowing for cosmological values for the Higgs boson mass in a window in 
which the electroweak vacuum is stable and therefore s ensitive to the field fluctuations during 
the early stages of the universe (jEspinosa et al. 20081 ). 



Limits of the validity of Higgs-type inflation have recently been debated by several 
authors. Specifically, Barbon & Espinosa (2009) argued that the large coupling of Higgs 
inflaton to the Ricci scalar makes this model invalid beyond the ultraviolet cutoff scale 
Ag ~ Mp/£ (here Mp = 2.4 x 10 18 GeV is the reduced Planck mass) which is below the Higgs 
field expectation value at N e- foldings during inflation, h ~ y/NMp/y/^. As consequence, at 
the ultraviolet cutoff scale A? at least one of the cross-sections of different scattering processes 
hits the unitarity bound ( iBurgess et al. 20091 ) . The fact that the quantum corrections due 
to the strong coupling to gravity makes the perturbative analysis to break down at energy 
scales above A^ was interpreted as a signature of a new physics, implying higher dimensional 
operators at energies above A^. However, the theory can still be considered valid above A^ 
if one finds some ultraviolet completion or if a very high degree of fine tuning is required, 
keeping in this way the unwanted contributions of higher dimensional operators small to zero 
( Bezrukov fc Shaposhnikov 2009[ De Simone et al. 2009 ). 



Recent papers (Lerner fc McDonald 2010al ; iLerner fc McDonald 2010bt IBurgess et al. 2010 



Hertzberg 20101 ) revisit the arguments against Higgs-type inflation addressing the issue of its 
naturalness with respect to perturbativity and unitarity violation in the Jordan and Einstein 
frames. It is shown that the apparent breakdown of this theory in the Jordan frame does 
not imply new physics, but a failure of the perturbation theory in the Jordan frame as a 
calculational method. These works demonstrate that for inflation based on a single scalar 
field with large non-minimal coupling, the quantum corrections at high energy scales are 
small, making the perturbative analysis valid. As consequence, for these models there is no 
breakdown of unitarity at the energy scale Ag. In particular, when the single- field Higgs 
inflation model is analyzed in the Einstein frame there is no breakdown of the theory at 
energy scales h > A^, where h is the canonically normalized Higgs scalar field in the Einstein 
frame. However, the inclusion of two or more scalar fields non-minimally coupled with gravity 
(in particular, the 3 Goldstone bosons of the Higgs d oublet) causes u nitarity violation in the 
Einstein frame at A^, making the theory unnatural (IHertzberg 2010l ). 



The present cosmological constraints on the Higgs mass are based on mapping be- 
tween the Renormalization Group (RG) flow equations and the spectral index of the cur- 
vature perturbations parameterized in terms of the number of e-foldings until the end 
of inflation, emerging from the analysis of CMB data combined with astrophysical dis- 
tance measurements. For a robus t interpretation of up coming observations from PLANCK 
( IMandolesi et al.2010f ) and LHC ( IBayatian et al. 20071 ) experiments it is imperative to un- 
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derstand how the inflationary dynamics of a non-minimally coupled Higgs scalar field may 
affect the degeneracy of the inflationary observables. 

The aim of this paper is to make a full proper analysis of the WM AP 7-year CMB mea- 



surements complem ented with astrophysical distance measurements (IKomatsu et al. 2010 



Larson et al. 20101 ) in the context of the non-minimally coupled Higgs inflaton field with 



gravity. The paper is structured as follows. In Section 2 we compute the power spectra of 
scalar and tensor density perturbations generated during inflation driven by a single scalar 
field non-minimally coupled to gravity. In Section 3 we derive the Higgs field equations 
and compute the RG improved Higgs field potential and in Section 4 we present our main 
results. In Section 5 we draw our conclusions. Throughout the paper a is the cosmological 
scale factor (ao = 1 today), k 2 = &ixM pl 2 where M p \ ~ 1.22 x 10 19 GeV is the present value 
of the Planck mass, overdots denotes the time derivatives and jip = d/d<p. 



2. COSMOLOGICAL PERTURBATIONS DRIVEN BY A 
NON-MINIMALLY COUPLED SCALAR FIELD 

In this section we compute the power spectra of scalar and tensor density perturbation 



generated during inf 
via the Ricci scalar ( 


ation driven by a single scalar field non-minimally coupled to gravity 


Fakir & Unruh 199(1 


Hwane & Noh 19961 


Komatsu & Futamase 1998; 


Komatsu & Futamase 1999; 


Hwane & Noh 2001: 


Tsuiikawa & Gumiudpai 20041) . The een- 


eral action for these models in the Jordan frame is given by f 


Futamase & Maeda 1989 





(1) 

where U((p) is a general coefficient of the Ricci scalar, R, giving rise to the non-minimal 
coupling, G((p) is the general coefficient of kinetic energy and V(ip) is the general potential. 
The generalized U(ip)H gravity theory in Equation (JT]) includes diverse cases of coupling. 
For the generally coupled scalar field U = (7 + K 2 ^ip 2 ), G(ip) = 1 and 7 and £ are constants. 
The non-minimally coupled scalar field is the case with 7 = 1 while the conformal coupled 
scalar field is the case with 7 = 1 and £ = 1/6. 

The conformal transformation for the action given in Equation ([T]) can be achieved by defining 
the Einstein frame metric as: 

9w = ng ll j,, £1 = 2k 2 U(<p), (2) 
where the quantities in the Einstein frame are marked by caret. 



Sj= d xy/^g 



%)R--%)(Vp) 2 -%) 
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The kinetic energy in the Einstein frame can be made canonical with respect to the 
new scalar field (p, defined through the scalar field propagator suppression factor s((p) as 
JPe Simone et al. 20091 ; iBarvinskv et al. 2009f ): 



1 + 3C/%), 



(3) 



2k 2 U 2 {ip) 

Thus the non-minimal coupling to the gravitational field introduces a modification to the 
Higgs field propagator by the factor s(<p), acting as back reaction of the gravitational field. 
The scalar potential V(ip) in the Einstein frame is given by: 

1 V(tp) 



leading to the following canonical form of the action in the Einstein frame: 



S, 



d 4 Xy/—g 



R (v^) 2 -y(£) 



2k 2 ' 2 



(4) 



(5) 



2.1. Background Field Equations 



When evaluating the field equations we assume that the background space-time can be 
written in the form of a flat (k=0) Robertson- Walker line element: 

ds 2 = g^ v dx^ dx u = -dt 2 + a 2 (t)dx 2 (6) 
= Q(x) (-di 2 + a 2 (t)dx 2 ) , 

where t is the cosmic time and a is the cosmological scale factor. From the above equation 
we obtain: 

da = Vttda, di=Vttdt. (7) 



Now the Friedmann equation in the Einstein frame can be written as (IKomatsu fc Futamase 1998 
Tsujikawa h Gumjudpai 2004 ): 



where: 



H 



2 _ k2 


[(f) 


3 




1 da 


1 


a dt 


y/Q 


dip 


(dtp 


dt 





H + 



1 dQ 
2Q~dt 



dt 



(8) 

(9) 
(10) 



Equations (Q and ffTUl) are enough to compute the background field evolution in the Einstein 
frame if the field equations in the Jordan frame are known (see the next section). 
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2.2. Scalar and Tensor Perturbations 



Neglecting the contribution of the decaying modes, the scale dependence of the ampli- 
tudes of scalar (S) a nd tensor (T) p erturbations in the Einstein frame are fully governed by 
the mode equation flMukanov 198ll ): 



d 2 u k 
dt 2 



+ k 



ld 2 z 
^~dP 



u k = 



where k is the comovi ng wave number of the mode function Uk- For the c ase of scalar 
perturbations we have (IHwang 19961 ; iHwang fc Noh 19961 ; iHwang fc Noh 20011 ) : 



1 d 2 z s 
z s dt 2 



[aHf 



(l + S s )(2 + 6 s + e) + -^r 

an 



'121 



where: 



zs = ay Qc 



Qs 



(13) 



and slow- roll parameters e and 5s are given by (jStewart Lyth 19931 ): 

e = -t, 5s = 4^-. (14) 
H 2 2HQ S 

In the case of tensor perturbations Equation (fl2l has the same form with the following 
replacements: 



z s -> z T = ay Q T , 



Qs Q T = 1 



0, 



T 



o 



(15) 



The power spectra of scalar and tensor perturbations are given by (ICopeland et al. 1994 ): 



V s (k) 



2?r 2 \Q 



i \ 2 I 19 

1 \ \u k \ 



V T {k) 



16fc 3 |-u fc | 



2 9 ' 

7rm£ a 2 



(16) 



and the spectral index of the scalar perturbations ris is obtained as usual as: n s — 1 
dlnPs(k)/dlnk. 



3. HIGGS BOSON AS INFLATON 



Higgs boson as inflaton adds non- minimal coupling to gravity (jBarvinskv &: Kamenshchik 1994: 
Bezrukov fc Shaposhnikov 2008 : Barvinsky et al. 2008 : De Simone et al. 20091 : Bezrukov et al. 2009 
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Fig. 1. — Left panel: The running of the coupling constants normalized to their initial values 
for m H = 132.8 GeV (continuous lines) with A(0) ~ 0.14 (blue), f (0) = 1.7 x 10 4 (magenta), 
m H = 160 (dashed lines) with A(0) = 0.21 (blue), f (0) = 2.1 x 10 4 (magenta) and m T = 171.3 
GeV with yt(0) = 0.91 (continuous red line). The green curves show the running of the Higgs 
field propagator suppression factor s(t). The right-hand gray region indicates the slow-roll 
inflationary regime. Right panel: The Einstein frame renormalization group improved po- 
tential as a function of the Higgs field for m# = 132.8 GeV (continuous line), m# = 160 
(dashed line) and rriTop = 171.3 GeV. In both cases we take the amplitude of scalar pertur- 
bations A s = 2.44 x 10~ 9 at the Hubble radius crossing £^=0.002 Mpc -1 and the vacuum 
expectation value v=246.22 GeV. 
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Bezrukov fc Shaposhnikov 20091 ). 

Taking the Higgs field potential V((p) of the Landau- Ginzburg type (by assuming that the 
spontaneous symmetry breaking arises through a condensate), the Jordan-frame effective 
action has the same form as given in Equation ([T]) with (see e.g. Futamase & Maeda 1989, 
Fakir & Unruh 1990, Makino & Sasaki 1991): 



U(<p) 



1 + k 2 & 2 

»-l/2 : 



V^) = U^-v 2 )'' 



G{ V ) 



(17) 



where: v = (v / 2GV)~ 1/2 =246.22 GeV is the vacuum expectation value of the Higgs field 
that sets the electroweak scale, A is the quadratic coupling constat of the Higgs boson with 
a mass m# = ^/2Xv and £ is the non-minim al coupling constant. The Jordan- frame field 
equations from the above action are given by (IKomatsu fc Futamase 19991 : iKaiser 19951 ): 



H 2 



K 



3(1 + K 2 i<\> 2 



V(<p) + -<p-6tH<p<p 



'11 



V + ^+{ 1 + K ^2 {1 + m )-- 1 + ^(1 + 60 ' (19) 

which in the slow- roll approximation ( \(p/(p\ <C H and |0 2 | V(<p) ) can be written as: 



K 



2 



H - w^m v(9) > (20) 

3 ^ " 1 + ^(1 + 60 • (21) 

The quantum corrections due to the interaction effects of the SM particles with Higgs boson 
through quantum loops modify the action coefficients U(<p), V((p) and G((p) from their 
classical expression given in Equations (1) an d (fTTl), taking the renormalization group (RG) 
improved forms UJt), V q (t), G q (t) defined as ( Barvinsky et al. 20081 : De Simone et al. 2009 



Clark et al. 2009: 



Lerner fc McDonald 2009h : 



U q (t) = ^- 2 (l + K%t)Gl{t)ip{t) 2 ) , (22) 

V q (t) = ^G 4 q (t)( V 2 (t)-v 2 (t)) 2 , (23) 
G q (t) = e-TW/a+TW) , (24) 

where j(t) is the Higgs field anomalous dimension given in the Appendix. 

The scaling variable t = ln^/m^) in the above equations normalizes the Higgs field and all 

the running couplings to the top quark mass scale m?- 
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As the energy scale of inflation is many order of magnitude above the electroweak scale 
(<p(t) » v), in the following we will approximate the Higgs potential by V((p) — A(/3 4 /4, 
neglecting the vacuum contribution and its running in the potential. 
Making the conformal transformation (J2J), Equations (j3J) and fll]) yield to: 

2 fd0(t)V i_ i + ^(ty 2 (t)(i + 6^)) 

{) W(t)J " 2k 2 (l + ^( t )^(t))2 ' 

v( f) J_ a (*)p 4 (*) , . 

1 J 16k 4 (1 + K 2 £(t)<p 2 (t)) 2 ' 1 ' 

The amplitude of scalar density perturbations at the Hubble radius crossing fc* is then given 
by: 



A 2 




(27) 



We compute the various t-dependent running constants, the Higgs field propagator suppres- 
sion factor and the Higgs field anomalous dimension by integrating the RG /3-functions as 
compiled in the Appendix. The runnings of SU{2) x 5(1) gauge couplings g',g, the SU(3) 
strong coupling g s , the top Yukawa coupling y t and the Higgs quadratic coupling A are com- 
puted by using two-loop quantum corrections while the running of non-minimal coupling 
constant £ is computed by using one-loop quantum corrections. 

One should note the importance of the quantum corrections due to non-minimal coupling. 
The quantum corrections to the classical kinetic sector G(ip) = 1 arise from the Higgs field 
anomalous dimension j(t) occurring with a factor of l/£ which in the inflationary regime 
(£ ~ 10 4 ) has a negligible small contribution. In the case of a classical gravity sector 
U(ip) = (1 + K 2 £(p 2 )/2K 2 , the conformal transformation (2) introduces a one-loop /3-function 
for £ with a term proportional to A due to Higgs running in a loop which has a small con- 
tribution during inflation due to the suppression of the Higgs fiel d propagator, while th e 



contribution of the remaining terms cancel to good approximation (IDe Simone et al. 20091 ). 
Although small, the one-loop quantum corrections due to the non-minimal coupling are not 
negligible but enough for the purpose of this analysis. 

For each case, the t-dependent running constants are obtained as: 

Y (t) = f t^tt dt', Y — {g, g>, g s , y t , A, £} , (28) 



t=o 1 + l{t') 

At t = 0, which corresponds to the top quark mass scale m T , the Higgs quadratic coupling 
A(0) and the top Yukawa coupling yt(0) are determined by the pole masses and the vacuum 
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expectation value v: 



A(0) = [1 + A H (m H )} , y t (0) = [1 + A T (m T )] , (29) 



where A H (m H ) and A T (m T ) are the corrections to Higgs and top quark mass respectively, 
computed following the scheme from the Appendi x of Espinosa et al. (2 008). 



The gauge coupling constants at uit scale are (jBarvinsky et al. 20091 ): g 2 (0) = 0.4202, 
g' 2 {0) = 0.1291 and g 2 (0) = 1.3460. The value of the non-minimal coupling constant £(0) is 
determined so that at the beginning of the slow-roll inflation tj n j the non-minimally coupling 
constant £(£j n i) is such that the calculated value of the amplitude of density perturbations 
given in Equation (j2"Tj) agrees with the measured value of A? s . 

Figure 1 presents the running of the coupling constants and of the Higgs field propagator 
suppression factor obtained for two different values of the Higgs boson mass. In both cases 
we also show the Einstein frame renormalization group improved potential as a function of 
the Higgs field ip = Ky/£<p(t). 



4. Results 



4.1. The CMB Angular Power Spectra 



We obtain the CMB temperature anisotropy and polarization power spectra by inte- 
grating the coupled Equations (jHD, ( |10|) and ( fill together with Equations ( 120]) and ( 12~TT) 
with respect to the conformal time imposing that the electroweak vacuum expectation value 
v =246.22 GeV is the true minimum of the Higgs potential at any energy scale (A(t) >0). 
We take wavenumbers in the range [5 x 10" 6 - 5] Mpc" 1 needed by the CAMB Boltzmann 
code (jLewis et al. 20001 ) to numerically derive the CMB angular power spectra and a Hubble 
radius crossing scale = 0.002Mpc~ 1 . The value of the Higgs scalar field at this scale 
is related to the quantum scale of infl ation <pr and to the du ration of inflation expressed in 
units of e- folding number N through (jBarvinsky et al. 20081 ): 



■ft 
<Pi 

Ai 



x 



NAi 

48^ ' 



^ (2s 4 + (g 2 + g' 2 ) 2 - IQyf) - 6A . 



(30) 
(31) 



where the inflationar y anomalous scaling parameter Ai (jBarvinsky fc Kamenshchik 1994 



Barvinsky et al. 20091 ) involves a special combination of quantum corrected coupling con- 
stants. These relations determine the value of the scaling parameter £* = ln^/mr) at 
Hubble radius crossing fc*. As the inflationary observables are evaluated at the epoch of 
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10 1 10 2 10 3 

multipole I 

Fig. 2. — The renormalization group improved CMB temperature and polarization angular 
power spectra (continuous lines) compared with the same power spectra obtained at the 
tree-level (dashed-lines) for m#=145 GeV (black lines) and m^=160 GeV (red lines). In 
both cases we take the top quark pole mass rnx = 171.3 GeV, the amplitude of scalar 
perturbations A 2 S = 2.44 x 10~ 9 at Hubble radius crossing ^=0.002 Mpc -1 and the vacuum 
expectation value v =246.22 GeV. 
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horizon-crossing quantified by the number of e-foldings N before the end of the inflation at 
which our present Hubble scale equalled the Hubble scale during inflation, the uncertainties 
in the determination of N translates into theoretical errors in determination of the infla- 



tionary observables (IKinney et al. 2004 ; Kinney &: Riotto 20061 ). Assuming that the ratio 
of the entropy per comoving interval today to that after reheating is negligible, the main 
uncertainty in the determination of N is given by the uncertainty in the determination of the 
reheating temperature after inflation. Recent studies of the reheating after inflation driven 
by SM Higgs field non-minimally coupled with gravity estimat es the reheating temperature 
in the range (jGarcia-Bellido et al. 20091 ; iBezrukov et al. 20091 ) : 



3.4 x 10 13 GeV < T r < 



A 



0.25 



1/4 



1.1 x 10 i4 GeV 



which translates into a negligible variation of the number of e-foldings with the Higgs mass 
(AN ~ 0.1). For the purpose of this work we choose N = k*/aH = 59 e-foldings in view 
of WMAP7+SN+BAO normalization at h ( komatsu et al. 20 id : [Larson et al. 201ol ). 
For each wavenumber k in the above range our code integrates the /3-functions of the t- 
dependent running constant couplings in the observational inflationary window imposing 
that k grows monotonically to the wavenumber k*, at the same time eliminating those mod- 
els violating the condition for inflation < e# = —H/H 2 < 1. 

Figure 2 presents the RG improved CMB temperature and polarization power spectra com- 
pared with the same power spectra obtained at the tree-level for m#=145 GeV and m#=160 
GeV. These plots clearly show that the CMB anisotropies are sensitive to the quantum ra- 
diative corrections of the SM coupling constants. 



4.2. Markov Chain Monte Carlo (MCMC) Analysis 



We use MCMC technique to reconstruct the Higgs field potential and to derive con- 
straints on the inflationar y observables and the Higgs mass from th e following datasets. 
The WMAP 7-year data (IKomatsu et al. 20101 ; lLarson et al. 20101 ) complemented with ge- 
ometric probes from the Type la supernovae (SN) distance-redshift relation and the baryon 
acoustic oscillations (BAO). The SN distance-redshift relation has been studied in detail 
in the recent u nified analysis of the published hetero geneous SN data sets - the Union 
Compilation08 fcowalski et al. 20081 : iRiess et al. 20091 ) . The BAO in the distribution of 
galaxies are extracted from Two Degr ee Field Galaxy Reds hidt Survey (2DFGRS) the Sloan 
Digital Sky Surveys Data Release 7 jPercival et al. 201ol ). The CMB, SN and BAO data 
(WMAP7+SN+BAO) are combined by multiplying the likelihoods. We use these mea- 
surements especially because we are testing models deviating from the standard Friedmann 
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Fig. 3. — The results of the fit of inflationary model with a non- minimal coupled Higgs 
scalar field to the WMAP7+SN+BAO dataset for top quark pole mass values: 168 GeV 
(green), 171.3 GeV (blue) and 173 GeV (magenta). The top plot in each column shows the 
probability distribution of different parameters while the other plots show their joint 68% 
and 95% confidence intervals. All parameters are computed at the Hubble crossing scale 
A;,=0.002Mpc- 1 . 
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expansion. These datasets properly enables us to account for any shift of the CMB angular 
diameter distance and of the expansion rate of the Universe. 

The likeli hood probabilities are evaluated by usin g: the public packages CosmoMC and 



CAMB ( ILewis Briddle 20021 ; iLewis et al. 20001 ) modified to include the formalism for 



inflation driven by non-minimally coupled Higgs scalar field as described in the previous 
sections. Our fiducial model is the ACDM standard cosmological model described by the 
following set of parameters receiving uniform priors: 

{fl b h 2 , Q c h 2 , 9 S , t, A 2 S , m H , m T ) , 

where: VL^h 2 is the physical baryon density, Q c h 2 is the physical dark matter density, S is the 
ratio of the sound horizon distance to the angular diameter distance, r is the reionization op- 
tical depth, A 2 S is the amplitude of scalar density perturbations, m H is the Higgs boson pole 
mass and ttlt is the top quark pole mass. For comparison we use the MCMC technique to 
reconstruct the standard inflation field potential and to derive constraints on the inflationary 
observables from the fit to WMAP7+SN+BAO dataset of the standard inflation model with 
minimally coupled scalar field. For this case we use the same set of input parameters with 
uniform priors as in the case of non-minimally coupled Higgs scalar field inflation, except 
for Higgs boson and top quark pole masses. The details of this computation can be found 
in Popa et al. (2009). 

For each inflation model we run 64 Mon te Carlo Markov chains , imposing for each case the 
Gelman & Rubin convergence criterion (jGelman fc Rubin 1992 ). Figure 3 presents the con- 



straints on the Higgs boson mass m#, the spectral index of the scalar density perturbations 
ns, the amplitude of the scalar density perturbations A 2 S and the ratio of tensor-to-scalar 
amplitudes R, as obtained from the fit to the WMAP7+SN+BAO dataset of the inflation 
model with non-minimally coupled Higgs scalar field for three different top quark pole mass 
values. We find that ns, A 2 S and R are dependent of the Standard Model parameters, in 
particular on the Higgs quadratic coupling and Yukawa coupling. One should recall that in 
the standard inflation these parameters are independent on the parameters of the Standard 
Model. 

The running of Higgs quadratic coupling A is increased for a heavier Higgs, also receiv- 
ing contributions from gauge couplings {g,g', g s } and top Yukawa coupling y t . In the in- 
flationary regime, the contribution from y t is increased as the top quark mass is varied 
toward higher mass values through its experimental allowed range: 168 GeV - 173 GeV 



( lAmsler et al. 20081 ). As a consequence, since we fixed the non- minimal coupling constant 



£ such that the amplitude of the scalar density perturbations A 2 S ~ A/£ 2 is at the observed 
value, A 2 S increases for a heavier Higgs boson and a higher top quark mass value, leading 
to the suppression of the spectral index of scalar density perturbations n s . Moreover, the 
joint confidence regions of the scalar spectral index n s and of the ratio of tensor-to-scalar 
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amplitudes R are anti-correlated. This can be attributed to a larger contribution of the ten- 
sor modes to the primordial density perturbations when Higgs boson and top quark masses 
are increased. Table 1 presents the mean values and the errors (68% CL) of the parameters 
from the posterior distributions obtained from the fit of the standard inflation model and 
the inflation model with non-minimally coupled Higgs scalar field with mT=171.3 GeV and 
^=246.22 GeV to WMAP7+SN+BAO dataset. We find for Higgs boson pole mass the fol- 
lowing dependence on m T and a s (m^)0 normalized in units of one standard deviations from 
their experimental central values: 



* (155.37 ± 3.85 ± 5) GeV + 3.8 GeV ( mT - 

V 2.3 Gov J 

- lAGev( as{mz) -° A176 ) (68% CL), (32) 



0.0020 

where we included the overall theoretical un certainty S ~ 2 GeV accounting for higher-order 



quantum corrections (jEspinosa et al. 20081 ) . In Figure 4 we present the dependence of the 
recovered Higgs field potential on m#, ris and R as obtained from the fit of inflationary 
model with non-minimally coupled Higgs scalar field to the WMAP7+SN+BAO dataset for 
different top quark pole mass values. Figure 4 explicitly demonstrates that the cosmological 
measurements not only probe the graviton-inflaton sector of the SM but also the variation 
of the scale of inflation due to the SM heavy particles coupled to inflation. 



CONCLUSIONS 



A number of papers have discussed bounds on the Higgs b oson mass coming fro m de- 
manding stability or metastability of the lifetime of the universe (jEspinosa et al. 20081 ). Fur- 
ther, by demanding that Higgs drive inflation, depending on the top quark mass and the com- 
putation of the RG improved effective potential, it was foun d that a heavier Higgs boson with 
a mass within the ab s olute stability bound s is required (IBezrukov fc Shaposhnikov 2 008; 
Barvinskv et al. 20081 : IBezrukov et al. 20091 : IBezrukov et al. 20091 : be Simone et al. 20091 : 
Bezrukov et al. 20091 ) . However, the present cosmological constraints on the Higgs boson 
mass are based on mapping between the RG flow and the scalar spectral index of of curva- 
ture perturbations. 



For a robu st interpretation of up coming observations from PLANCK (IMandolesi et al.2010f ) 
and LHC (IBayatian et al. 20071 ) experiments it is imperative to understand how the infla- 
tionary dynamics of a non-minimally coupled Higgs scalar field with gravity may affect the 



l a s = gg/4:ir is the effective QCD coupiing constant 
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Fig. 4. — The dependence of the reconstructed Higgs field potential (the joint 68% and 
95% confidence intervals) on m#, n$ and R as obtained from the fit of inflationary model 
with non-minimally coupled Higgs scalar field to the WMAP7+SN+BAO dataset for top 
quark pole mass values: 168 GeV (green), 171.3 GeV (blue) and 173 GeV (magenta). All 
parameters are computed at the Hubble crossing scale /c^=0.002Mpc _1 . 
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determination of the inflationary observables. The aim of this paper is to make a full proper 



analysis of the WMAP 7-year CMB measurements ( IKomatsu et al. 2010l ; lLarson et al. 20101 ) 



complem ented with geometric probes from the Type la supernovae (SN) distance- redshift 



relation (IKowalski et al. 20081 ; iRiess et al. 20091 ) and the baryon acoustic oscillations (BAO) 



in the distribution of galaxies from Two Degree Fie ld Galaxy Redshidt Survey (2DFGRS) 



and the Sloan Digital Sky Surveys Data Release 7 (jPercival et al. 20101 ). in the context of 
the non-minimally coupled Higgs inflaton with gravity. 

We compute the full RG improved effective potential including two-loop beta functions for 
SU(2) x £7(1) gauge couplings g',g, the S77(3) strong coupling g s , the top Yukawa coupling 
y t and the Higgs quadratic coupling A and one-loop beta functions for non-minimal coupling 
constant £ and vacuum expectation value v. We also include the curvature in RG flow equa- 
tions through Higgs field propagator suppression function s(t) and the Higgs field anomalous 
dimension ^(t). 

The initial conditions for A and yt are properly obtained through the pole mass matching 
scheme while the inflationary anomalous scale parameter Ai relates the initial value of the 
Higgs inflation field to the quantum scale of inflation and the number of e-foldings. 
We use MCMC technique to reconstruct the Higgs field potential and to derive constraints 
on the inflationary observables and the Higgs mass from WMAP7+SN+BAO dataset. For 
the central value of the top quark pole mass my = 171.3GeV the fit of the inflation model 
with non-minimally coupled Higgs scalar field to WMAP7+SN+BAO dataset leads to the 
following 95% CL bounds on Higgs boson mass: 

143.7 GeV < m H < 167.0 GeV , (33) 

where we take into account the overall theoretical error 5 = ±2 GeV and a s (mz) = 0.1176. 
We show that the inflation driven by a non-minimally coupled scalar field to the Einstein 
gravity leads to significant constraints on the scalar spectral index n$ and tensor-to-scalar 
ratio R, when compared with the similar constraints from the standard inflation with min- 
imally coupled scalar field. In particular, one should note the smallness of tensor-to-scalar 
ratio (R ~ 10 -3 ) that is challenging the future polarization experiments. 
We conclude that in order to obtain an accurate reconstruction of the Higgs potential in 
terms of inflationary observables it is imperative to improve the accuracy of other parame- 
ters of the SM as the top quark mass and the effective QCD coupling constant. 
For example, it is expected that in the near future LHC will improve the determination of 
the current value of top quark mass to Am T ~ 1.5 GeV. From Equation ( I3"2"j) it follows that 
this improvement will lead to an improvement in the determination of the Higgs boson mass 
to AmT — 2.4 GeV. Since A 2 S ~ A, using Equation (12?]) with R = 16e and fixing all param- 
eters at their observed values, it follows that the expected improved determination of the 
top quark mass leads to an improved accuracy in the determination of the Higgs potential 
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of about 3%. 
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6. APPENDIX 



In this appendix we collect the SM renormalization group /?- functions (IFord et al. 19921 ). 
including the Higgs field propagator suppression factor s(t) given in Equation ( |25|) . at the 
renormalization energy scale t = ln(ip/mt) beyond the top quark mass mt. 



The two-loop /3-functions for gauge couplings Qi = {g',g,g s } are ( lEspinosa et al. 20081 ): 



0: 



2^3 



kgfbi + k l g 



J2 B v9* ~ s(t)<%y 2 t 
.i=i 



(34) 



where k = 1/167T 2 and 



b = ((40 + s(t))/6, -(20 - s(t))/6, -7), B 

d l = (17/6,3/2,2). 



199/18 9/2 44/3 
3/2 35/6 12 
11/6 9/2 -26 



(35) 



For the top Yukawa coupling y t , the two- loop /3-f unction is given by (IDe Simone et al. 20091 ): 



tit 



ky t 



+ k 2 y t 



23 4 3 2 /2 i H87 /4 2 2 

- T -~ A 9 9 +^g +$gg s 



+ fg' 2 g 2 s - W8gt + {^g 2 + ^g' 2 + 36g 2 ) s{t)y 2 t + 6 (-2s 2 (t)yt - 2s\t)y 2 \ + S 2 (t)A 2 ) 



(36) 



The two-loop ft- function for the Higgs quadratic coupling A is ( IDe Simone et al. 20091 ): 

3 



Px = k 
+ k 2 



24s 2 A 2 - 6yf + | (V + (<? 2 + </ 2 ) 2 ) + (-9<? 2 - 3</ 2 + 12y 2 ) A 



1 (915/ - 289^ - 559sV 4 - 379</ 6 ) + 30s(t)y« - yf + 32g 2 + 3s(t)A 



,/2 
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Table 1: The mean values from the posterior distributions of the parameters obtained from 
the fit of the standard inflation model and Higgs inflation model with mr=171.3 GeV and 
v=246.22 GeV to WMAP7+SN+BAO dataset. The errors are quoted at 68% CL. All 
parameters are computed at the Hubble radius crossing A;*=0.002 Mpc" 1 . 



Model 


Standard Inflation 


Higgs Inflation 


Parameter 






ioon b h 2 


2.259±0.054 


2.257±0.051 


n c h 2 


0.113±0.003 


0.114±0.003 


T 


0.088±0.015 


0.086±0.013 


S 


1.038±0.002 


1.037±0.002 


ln[10 10 A|] 


3.157±0.031 


3.161±0.032 


n s 


0.960±0.012 


0.972±0.0004 


R 


< 0.144 


0.0036±0.0009 


m H (GeV) 




155.372±3.851 


X 




0.216±0.053 


£ x io- 4 




3.147 ± 0.509 
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A (-f g 4 + ^g 2 g' 2 + f^(*)</ 4 + 108*WA + 36s 2 (ty 2 A - 312 s 4 (t)A 2 ) 



9 



21 



19 



Vt ( —X + -9 2 9 12 ~ ^ + A ( -7^9 + Y 9' 2 + 80s 2 - 144s 2 (t)A 



45 



85 



(37) 



The one-loop /^-function for non-minimal cou pling £ is given by (IBezrukov fc Shaposhnikov 2009 ; 



Clark et al. 20091 : iLerner fc McDonald 20091 ): 

* = k (* + l) ( 6(1 + S2(m + 6V * ~ I 9 '' " I 9 ') ■ 



(3? 



The reference Bezrukov & Shaposhnikov (2009) also gives the one-loop /3-function for the 
vacuum expectation value v 2 in the form: 



A, 2 = M^ ,2 + 35 2 -6^K- 



(39) 



Finally, the two- loop Higgs field anomalous dimension 7 is given by ( )De Simone et al. 20091 ) : 

'9g 2 , 3g' 2 



7 = —k 

+ * 2 



4 + 4 



45 



-3^ 



85 



- k 



271 4 __9_ 2 , 2 _431 
32 ^ 16 P 96 S[ )9 



27 



/■i 



g 2 + -g' 1 + 20<? 2 ) yf + —s(t)tf - 6s d (t)A 



(40) 
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